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Abstract. We utilize a sample of galaxy clusters at 0.35<,z<0.6 drawn from the Las Cam- 
panas Distant Cluster Survey (LCDCS) to provide the first non-local constraint on the cluster- 
cluster spatial correlation function. The LCDCS catalog, which covers an effective area of 
69 square degrees, contains over 1000 cluster candidates. Estimates of the redshift and velocity 
dispersion exist for all candidates, which enables construction of statistically completed, volume- 
limited subsamples. In this analysis we measure the angular correlation function for four such 
subsamples at z~0.5. After correcting for contamination, we then derive spatial correlation 
lengths via Limber inversion. We find that the resulting correlation lengths depend upon mass, 
as parameterized by the mean cluster separation, in a manner that is consistent with both local 
observations and CDM predictions for the clustering strength at z=0.5. 



1 Introduction 

The spatial correlation function of galaxy clusters provides an important cosmological test, as both 
the amplitude of the correlation function and its dependence upon mean intercluster separation are 
determined by the underlying cosmological model. In hierarchical models of structure formation, the 
spatial correlation length, ro, is predicted to be an increasing function of cluster mass, with the precise 
value of ro and its mass dependence determined by <7g (or equivalently Qo> using the constraint on 
o"8 — r?o from the local cluster mass function) and the shape parameter T. Low density and low 
r models generally predict stronger clustering for a given mass and a greater dependence of the 
correlation length upon cluster mass. 

In this paper we utilize the Las Campanas Distant Cluster Survey (LCDCS) to provide a new, 
independent measurement of the dependence of the cluster correlation length upon the mean inter- 
cluster separation (d c = n 1 / 3 ) at mean separations comparable to existing Abell and APM studies. 
We first measure the angular correlation function for a series of subsamples at z ~ 0.5 and then derive 
the corresponding ro values via the cosmological Limber inversion |], 10, The resulting values 



constitute the first measurements of the spatial correlation length for clusters at z ^ 0.2. Popular 
structure formation models predict only a small amount of evolution from z = 0.5 to the present - a 
prediction that we test by comparison of our results with local observations. 



2 The Survey 

The recently completed Las Campanas Distant Cluster Survey is the largest published catalog of 
galaxy clusters at z ^ 0.3, containing 1073 candidates ||. Clusters are detected in the LCDCS as 
regions of excess surface brightness relative to the mean sky level, a technique that permits wide-area 
coverage with a minimal investment of telescope time. The final statistical catalog covers an effective 
area of 69 square degrees within a 78° x 1.6° strip of the southern sky (860 x 24.5 h^ 1 Mpc at z=0.5 
for $7o=0.3 ACDM). Gonzalez et al. (2001a) also provide estimated redshifts (z es t), based upon the 
brightest cluster galaxy (BCG) magnitude-redshift relation, that are accurate to ~15% at z es t = 0.5, 
and demonstrate the existence of a correlation between the peak surface brightness, S, and velocity 
dispersion, a. Together these two properties enable construction of well-defined subsamples that can 
be compared directly with simulations and observations of the local universe. 
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Figure 1: a/ Extinction-corrected surface brightness distribution of the LCDCS clusters used to 
construct subsamples. Gonzalez et al (20001a) found that logo" oc logE(l + z) 5 .l, and so we employ 
this relation to definty subsamples that are approximatlye dispersion-limited. The solid lines denote 
the redshift and E(l + z) 5 ' 1 thresholds of the subsamples. The dotted curve is the surface brightness 
threshold of the full LCDCS catalog for F>(B — V)=0.05. Errors in the estimated redshifts move 
individual data points along tracks similar to this curve; b/ Angular correlation functions for the 
entire LCDCS catalog and for a volume-limited subsample with 0.35 < z < 0.475. The solid lines are 
best-fit power law models. 



2.1 The LCDCS Angular Correlation Function 

To compute the two-point angular correlation function, we use the estimator of Landy & Szalay (1993). 
We measure the angular correlation function both for the full LCDCS catalog and for three approx- 
imately velocity dispersion-limited subsamples at z~0.5 (Figure |l]a). We restrict all subsamples to 
z es t > 0.35 to avoid incompleteness, while the maximum redshift is determined by the surface brightness 
threshold of the subsample. 

The angular correlation function for the entire LCDCS catalog is shown in the upper panel of 
Figure [l]&, with logarithmic angular bins of width 5 log 6=0.2. Modeling this correlation function as 
a power law, ui{6) = A^O 1 ' 1 = (#/6>n) 1-7 , a least-squares fit for all LCDCS clusters over the range 
2'-5° yields 7=1.83±0.12 and #n=56±22". The angular correlation function for the lowest redshift 
subsample is shown in the lower panel of Figure [Tjfcf, overlaid with a best-fit power law. We derive 
best-fit values both allowing 7 to vary as a free parameter and also fixing 7=2.1 — equivalent to the 
best-fit value for the lowest redshift subsample and the best fit value for the ROSAT All-Sky Survey 



1 Bright Sample [13]. 



We then apply a correction to these best-fit values to account for the impact of false detections in 
the LCDCS catalog, which for this data set act to suppress the amplitude of the observed correlation 
function. If we assume that the contamination is spatially correlated and can be described by a 
power law with the same slope as the cluster angular correlation function (a reasonable approximation 
because for galaxies - which are likely the primary contaminant - 7~1.8-1.9 [e.g. 2, 19]), then the 
observed value of A w is 

= A:luster(l ~ f) 2 + ^.false/ 2 , (1) 

where / is the fractional contamination. For detections induced by isolated galaxies of the same 
magnitude as BCG's at z ~ 0.35 (and identified as galaxies by the automated identification criteria 
described in Gonzalez et al. (2001a), we measure that A ga \ is comparable to A LJA , the net clustering 
amplitude for all LCDCS candidates at 0.3<z es t<0.8. For detections identified as low surface bright- 
ness galaxies (including some nearby dwarf galaxies) we measure ^4lsb ~ IOA^a • While these systems 
are strongly clustered, we expect that they comprise less than half of the contamination in the LCDCS. 
For multiple sources of contamination the effective clustering amplitude ^4f a i S e = /i) 2 > so 



Table 1: Spatial Correlation Lengths 



~z ACDM (fi = 0.3) OCDM (fi = 0.3) TCDM(r=0.2) 

range d c r d c r d c r 



0.35-0.475 38.4 15.1±|;g 33.8 13.2±^;| 30.9 12.1; 
0.35-0.525 46.3 18.5±|;| 40.6 16.1±|;| 36.9 14.7; 
0.35-0.575 58.1 22.1 j|;g 50.8 19.1±| : | 46.0 17.3; 



Note — Units of ro and d c are h 1 Mpc. 



the effective clustering strength of the contamination is Af a i sc <^ 2.5A ulA even including the LSB's. 



3 The Spatial Correlation Length 

The observed angular correlation function can be used to determine the three-space correlation length 
if the redshift distribution of the sample is known. This is accomplished via the cosmological Limber 
inversion || [l(|. For a power-law correlation function with redshift dependence f(z), 



f(r) = - 



The corresponding comoving spatial correlation length is ro(z) = rof(z) 1 /~ / , and the Limber equation 



is 



'I = A, 



r( 7 /2) 



fr r(i/2)r[( 7 -i)/2] 



J z f(dN/dz) 2 E(z)D A (z) 1 ^f(z)(l + z)dz 



f£(dN/dz)dz 



(3) 



where dN/dz is the redshift distribution of the sample, Da(z, f^cb ^a) is the angular diameter distance, 
and E(z) is defined as in Peebles (1993). Because little evolution in the clustering is expected over the 
redshift intervals spanned by our subsamples (see Figure ||), f(z) can be pulled out of the integral. 

For the LCDCS we estimate the true redshift distributions of our subsamples based upon the 
observed distribution of estimated redshifts. If we approximate the redshift error distribution as 
Gaussian with a z /z ~ 0.14 at z=0.5 ||, then the actual redshift distribution dN/dz is approximately 
equal to the observed redshift distribution dN b s /dz for a given subsample convolved with this Gaussian 
scatter. To test the validity of this approach, we also try modeling dN/dz using the theoretical mass 
function of Sheth & Tormen (1999) convolved with redshift uncertainty. Comparing these two methods 
we find that the derived spatial correlation lengths agree to better than 3% for all subsamples .[] 



4 Results and Comparison with Local Data 

Table |l| lists the correlation lengths (ro) and mean intercluster separations (d c ) that we derive for 
the three LCDCS subsamples. The values of ro and d c are cosmology-dependent, so we list these 
quantities for three different cosmologies — ACDM (O =0.3,r=0.2), OCDM (O =0.3,r=0.2), and 
rCDM (r=0.2). We opt to fix 7=2.1 when deriving the correlation lengths in Table [l] because 7 is 
strongly covariant with A u and hence poorly constrained by the LCDCS data set. In no instance does 
this choice alter the derived ro value by more than 10% from the value obtained when 7 is treated as 
a free parameter, but we caution that it can systematically bias the observed dependence of ro upon 
d c . For subsamples with best-fit values of 7>2.1, fixing gamma slightly increases the derived value of 
ro, which results in a mild steepening of the dependence of ro upon d c for the LCDCS subsamples. 



This result is insensitive to the exact mass threshold assumed for the theoretical mass function. 
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Figure 2: a/ Comparison of the LCDCS data with local samples and theoretical predictions for 
ACDM (f)=0.3). The error bars on the LCDCS data correspond to the 1-a statistical uncertainty. 
Overlaid lines correspond to analytic predictions (T=0.2) for z=0 and z=0.5, and to results from the 
Virgo Consortium Hubble Volume simulations; b/ Same plot for rCDM. 



The cluster correlation function is predicted to not evolve significantly between z=0.5 and the 
present; we thus compare the LCDCS results directly with local observations, as is shown in Figure ^ 
for ACDM and rCDM. For both cosmologies the LCDCS values of rn are comparable to those from 
the Edinburgh-Durham Galaxy Catalogue [14], APM survey EJ, and the MX Survey northern sample 
[O], but smaller than those found by Peacock &; West (1992) for the Abell catalog. Also shown are 
the lowest d c data points for the XBAC catalog [ffl, which probe higher masses than our study. 

We also plot theoretical predictions for comparison with the observational data. Results from the 
Virgo Consortium Hubble Volume simulations are shown as a dot-dash line in Figure |2]a || . The other 
lines are analytic predictions based upon the work of Sheth & Tormen (1999). Like the APM data, 
the LCDCS results are consistent with the low-density models (independent of £Ia)- Li contrast, the 
plotted rCDM model systematically underestimates both the local and z=0.5 data. rCDM can be 
made to match the data only by decreasing T to values that are inconsistent with constraints from 
the galaxy power spectrum Q. 

To assess the robustness of these results, we also quantify potential systematic biases (see Gonzalez 
et al. 20016 for more detail). We find that our results can only be significantly altered if the uncertainty 
in the estimated redshifts is underestimated. If so, then the correlation lengths we derive would be 
systematically too small (by ~2 h~ l Mpc if a z j z=0.25, for example). The second most significant 
potential systematic arises from fixing 7, which as mentioned above may lead us to overestimate the 
dependence of r$ upon d c . Treating 7 as a free parameter would reduce the derived r$ values for the 
three subsamples (in order of increasing d c ) by 0%, 7%, and 10%, but would not qualitatively change 
our results. 



5 Discussion and Conclusions 

The Las Campanas Distant Cluster Survey is the largest existing catalog of clusters at z>0.3, providing 
a unique sample with which to study the properties of the cluster population. We use the LCDCS 
to constrain the cluster-cluster angular correlation function, providing the first measurements for a 
sample with a mean redshift z^0.2. From the observed angular correlation function, we derive the 
spatial correlation length, rn, as a function of mean separation, d c . Only modest evolution in the 
clustering amplitude is predicted between z = 0.5 and present, and we thus compare our results 
directly with local data. We find that the LCDCS correlation lengths agree with results from local 
samples, and observe a dependence of rn upon d c that is comparable to the results of Croft et al. 



(1997) for the APM catalog. This clustering strength, its dependence on number density, and its 
minimal redshift evolution are consistent with analytic expectations for low density models, and with 
results from the ACDM Hubble Volume simulations. Consequently, while statistical uncertainty limits 
our ability to discriminate between cosmological models, our results are in concordance with the flat 
ACDM model favored by recent supernovae and cosmic microwave background observations H IHi HO] • 
A final result of this analysis is that it demonstrates the utility of large catalogs like the LCDCS 
that are statistical in nature. While the properties of any particular cluster in the LCDCS catalog 
are rather uncertain, the properties of the sample as a whole are well-defined, which is sufficient 
for constraining properties such as the clustering strength and evolution in the comoving number 
density. This statistical approach requires a relatively small investment in telescope time, and so can 
be extended in the future to much larger samples than the LCDCS. 
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